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ABSTRACT 

Recent studies have shown that the extended main-sequence turn off (eMSTO) is a common 
feature of intermediate-age star clusters in the Magellanic Clouds (MCs). The most simple 
explanation is that these stellar systems harbor multiple generations of stars with an age dif¬ 
ference of a few hundred Myrs. 

However, while an eMSTO has been detected in a large number of clusters with ages 
between ~l-2 Gyrs, several studies of young clusters in both MCs and in nearby galaxies 
do not find any evidence for a prolonged star-formation history, i. e. for multiple stellar gen¬ 
erations. These results have suggested alternative interpretation of the eMSTOs observed in 
intermediate-age star clusters. The eMSTO could be due to stellar rotation mimicking an age 
spread or to interacting binaries. In these scenarios, intermediate-age MC clusters would be 
simple stellar populations, in close analogy with younger clusters. 

Here we provide the first evidence for an eMSTO in a young stellar cluster. We ex¬ 
ploit multi-band Hubble Space Telescope photometry to study the ~300-Myr old star cluster 
NGC 1856 in the Large Magellanic Cloud and detected a broadened MSTO that is consistent 
with a prolonged star-formation which had a duration of about 150 Myrs. 

Below the turn-off, the MS of NGC 1856 is split into a red and blue component, hosting 
33+5% and 67+5% of the total number of MS stars, respectively. We discuss these findings in 
the context of multiple-stellar-generation, stellar-rotation, and interacting-binary hypotheses. 

Key words: globular clusters: individual (NGC 1856) — stars: Population II 


1 INTRODUCTION 

The finding that intermediate-age star clusters in the Magellanic 
Clouds exhibit a bimodal or extended main-sequence turn off (eM¬ 
STO) is one of the most intriguing discoveries in the field of stellar 
populations of the last decade (e.g. Bertelli et al. 2003; Baume et 
al. 2007; Mackey & Broby Nielsen 2007; Glatt et al. 2008; Goud- 
frooij et al. 2011; Keller et al. 2012; Rubele et al. 2013). High- 
accuracy photometry with the Hubble Space Telescope (HST) has 
revealed that the eMSTO is a common feature among the ~l-2 Gyr- 
old MCs’ clusters as it has been detected in most of the objects stud¬ 
ied so far (Milone et al. 2009 — hereafter Paper I —; Goudfrooij et 
al.2014). 


The most straightforward interpretation of the eMSTO is that 
these stellar systems have experienced an extended star-formation 
history with a duration of 100-500 Myr (e.g. Mackey et al. 2008; 
Paper I; Goudfrooij et al. 2009, 2014). As an alternative explana¬ 
tion, Bastian & de Mink (2009) and Yang et al. (2011, 2013) have 
suggested that stellar rotation or interacting binaries can mimic an 
age spread and could be responsible for the eMSTO (but see Girardi 
et al. 2011 and Platais et al. 2012). 

Noticeably, there is no evidence for prolonged star formation 
in clusters younger than ~1 Gyr although age spreads have been 
searched in a large number of young massive clusters with differ¬ 
ent techniques (Larsen et al. 2011; Milone et al. 2013 — hereafter 
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Paper II —; Bastian & Silva-VilIa20I3; Bastian et al. 2013; Nieder- 
hofer et al. 2015). 

The properties of the progenitors of the clusters where we ob¬ 
serve eMSTO have been widely studied by Bekki & Mackey (2009) 
and Keller et al. (2011) in the context of a prolonged star-formation 
history. A sufficiently large mass should enable a cluster, in prin¬ 
ciple, to retain ejecta from a previous stellar generation which pol¬ 
luted the gas from which the second generation can form. This pos¬ 
sibility implies that the mass of the progenitor should have been 
of the order of 10^Al© at the time of cluster formation. Keller 
et al. (2011) concluded that NGC 1856 in the Large Magellanic 
Cloud, with a mass of ~ 10^ AT© (Mackey & Gilmore 2003) and an 
age of ~300 Myr (Bastian & Silva-Villa2013) is the only cluster 
with measured age and mass close to that of an eMSTO progenitor. 

NGC 1856 has been previously studied by Bastian & Silva- 
Villa(2013) who have analyzed Wide-Field-Planetary Camera 2 
(WFPC2) photometry from Brocato et al. (2001) and found an age 
of 280 Myr and no evidence for an eMSTO and an age spread. 
In this paper we exploit a new dataset, which consists of multi¬ 
wavelength images collected with the Ultraviolet and Visual chan¬ 
nel of the Wide Field Camera 3 (UVIS/WFC3) on board of HST, 
to further investigate NGC 1856 in the context of multiple stellar 
populations. 


2 DATA AND DATA ANALYSIS 

The dataset used in this paper is summarized in Table [T] and con¬ 
sists of images collected with UVIS/WFC3 through five filters that 
cover a wavelength range from the ultraviolet to ~8100 A. The poor 
charge-transfer efficiency (CTE) in the UVIS/WFC3 images have 
been corrected by using a software written by Jay Anderson fol¬ 
lowing the recipe by Anderson & Bedin (2010). 

Photometry and astrometry of the images have been car¬ 
ried out with a software program mostly based on the Anderson 
et al. (2008) software adapted to UVIS/WFC3 by Jay Anderson. 
Briefly, we used two distinct methods to infer luminosity and posi¬ 
tion of bright and faint stars. Astrometry and photometry of bright 
stars have been performed independently in each image by us¬ 
ing the best point-spread function (PSF) model available, and later 
combined. For that purpose, we used library PSF by Anderson et 
al. (in preparation) and accounted for small focus variations due to 
the ‘breathing’ of HSThy deriving a spatially-constant perturbation 
for each exposure. The flux and the position of each very faint stars 
can be determined more robustly by simultaneously fitting all the 
pixels in all the exposures. Moreover, we adopted the solution of 
Bellini & Bedin (2009) and Bellini, Anderson & Bedin (2011) to 
correct stellar positions for geometrical distortion. The photometry 
has been calibrated as in Bedin et al. (2005) and by using the zero 
points provided by the STScI web page for WFC3/UVI^ 

Since the study of multiple populations in star clusters requires 
high-precision photometry, we have used a number of indexes pro¬ 
vided by the adopted software as diagnostics of photometry quality 
(see Anderson et al. 2008 for details). Similarly with what was done 
in Papers I (see their Sect. 2.1) and II, we have limited our study of 
NGC 1856 to a sub-sample of stars that have small astrometric er¬ 
rors, are relatively isolated, and well fitted by the PSF. 

We show in Fig.[T]the trichromatic stacked images of the ana¬ 
lyzed field of view (left panel) and the mpssew vs. raptsew “ ™f 814 w 

* http://www.stsci.edu/hst/acs/analysis/zeropoints/zpt.py 


CMD of all the stars for which photometry in both F336W and 
F814W bands is available (right panel). The two circles superim¬ 
posed on the left-panel image have a radius of 25 arcsec and delimit 
a region named cluster field which is centered on NGC 1856 and is 
mostly populated by cluster members (green circle), and a refer¬ 
ence field, where the contamination from NGC 1856 stars is negli¬ 
gible (yellow circle). Stars in these two regions will be used in the 
study of the multiple stellar populations of NGC 1856 in Sects. HIS] 

2.1 Artificial Stars 

We have used the method described in details by Anderson et 
al. (2008) to perform the artificial-star (AS) tests. Briefly, a list of 
10^ stars has been generated and placed along the fiducial line of 
the MS of NGC 1856. The list includes the coordinates of the stars 
in the reference frame and the magnitudes in all the bands studied 
in this paper. ASs have been distributed across the field of view 
according to the overall cluster distribution as in Paper 1. 

In each image and for each star in the input list, the software 
by Anderson et al. (2008) generates a star with appropriate flux and 
position and measures it by using the same procedure as for real 
stars. 

For the ASs the software provides the same diagnostics of the 
photometric quality as for real star. We applied to ASs the same 
procedure as used for real stars to select a sub-sample of relatively- 
isolated stars with small astrometric errors, and well fitted by the 
PSF. ASs have been used to estimate errors of the photometry and 
estimate the fraction of chance-superposition binaries. Moreover, 
ASs have been used to derive the completeness level of our sample 
as in Paper 1. 


3 THE CMD OF NGC 1856 

The rap 336 w vs. mF 336 w “ ™f 814 w CMD of stars in the cluster field 
is plotted in the upper-left panel of Fig. It shows that the MS 
and the MSTO of NGC 1856 are widely spread in color and mag¬ 
nitude. A visual comparison between the observed CMD and the 
error bars, derived from AS tests and plotted on the left, reveals 
that this spread can not entirely be explained by the photometric 
errors. In the following we investigate whether this broadening is 
intrinsic or can be due to differential reddening, spatial variation of 
the photometric zero point across the field of view, field stars, or 
binaries. 

3.1 Differential reddening 

In order to investigate the effect of differential reddening on the 
MSTO of NGC 1856, in the upper-right panels of Fig. |2]we have 
plotted the CMD not corrected for differential reddening for the 
stars located in six circular sub-regions (I-VI) of 25-arcsec radius. 
Due to the small area of each region, the variation of reddening 
therein is smaller than that of the entire cluster field. Noticeably 
the color and magnitude of the MSTO is broadened in each CMD, 
thus suggesting that the young NGC 1856 hosts an eMSTO in close 
analogy with intermediate-age MC clusters. 

To properly address this issue, we have corrected the photome¬ 
try of stars in the cluster field for differential reddening as in Milone 
et al. (2012). Briefly, we started to select a sample of reference stars 
including those bright MS stars which pass the criteria of selection 
described in Sect. |2] Then we have determined the fiducial line of 
these reference stars and calculated, for each of them, the distance 
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Figure 1. Left panel: Stacked trichromatic image of the field around NGC 1856 studied in this paper. Green and yellow circles mark the cluster and the 
reference field, respectively. Right panel: Hip 33 gw vs.mp 336 w “ «tF8l4W CMD of all the stars in the analyzed field of view for which F336W and F814W 
photometry is available. 


Table 1. Description of the data set used in this paper. 


INSTR. 

DATE 

NxEXPTIME 

FILTER 

PROGRAM 

PI 

WFC3AJVIS 

Feb 09, Mar 24, May 18, and Jun 11, 2014 

11x711s 

F336W 

13379 

A. P. Milone 

WFC3/UVIS 

Nov 12 2013 

185s+2x430s 

F438W 

13011 

T. H. Puzia 

WFC3/UVIS 

Nov 12 2013 

2x350s 

F555W 

13011 

T. H. Puzia 

WFC3/UVIS 

Nov 12 2013 

735s+2x940s 

F656N 

13011 

T. H. Puzia 

WFC3/UVIS 

Nov 12 2013 

51s+360s+450s 

F814W 

13011 

T. H. Puzia 

WFC3/UVIS 

Feb 09, Mar 24, May 18, and Jun 11, 2014 

4x90s+4x704s 

F814W 

13379 

A. P. Milone 


from the fiducial line along the reddening line. The reddening di¬ 
rection has been determined by using the absorption coefficient for 
UVIS/WFC3 filters kindly provided by Aaron Dotter (2014, pri¬ 
vate communication). Namely these are jE(B - V) = 5.10, 

Af 438 w/£(B -V) = 4.18, AF555w/£(fi - V) = -i-ll, Af 656 n/£(S - 
V) - 2.54, and Af 814 w/F(B - V) = 1.86. As differential reddening 
of each star in the cluster field we assumed the median distance of 
the 45 closest reference stars. The uncertainty on the differential 
reddening has been estimated as the ratio between the dispersion of 
the 45 values of the distance of each reference star from the fidu¬ 
cial line and the square root of 44. We have found that the redden¬ 
ing in the direction of NGC 1856 is not uniform and ranges from 
AE{B - V) —0.05 to —hO.07 with respect to the average reddening 
E(B - V)=0.15 (see Sect.lSH. 

The mF335w vs. raF 336 w-™F 8 i 4 w CMD corrected for differential 
reddening is shown in the lower-left panel of Fig. [2] The error bars 
plotted on the left account also for the uncertainty in the differential 
reddening correction above estimated. The corrected photometry 
unveils several features in the CMD of NGC 1856. 

The color width of the MS, below /nF 336 w ~21.5, is compara¬ 
ble to the measurement errors and a large number of binary systems 
is visible on the red side of the MS. The MSTO is widely spread in 


color and magnitude thus indicating that NGC 1856 belongs to the 
class of clusters with an eMSTO. Furthermore, we confirm previous 
findings by Bastian et al. (2013) that the red clump of NGC 1856 
spans a wide color range and hosts two main components clustered 
around (mF 336 w - mF8i4w:«rF336w)=(2.60:19.25) and (2.75:19.25). 

For completeness, we show the reddening map in the lower- 
right panel of Fig.|2l where we also mark the sub-regions I-VI, used 
to derive the six CMDs in the upper-right panel. 

In order to investigate the MSTO morphology of NGC 1856 
we have used the mF 335 w vs. wtFtsew “ ™f 555 w CMD for stars in 
the cluster field of view corrected for differential reddening (panel 
(ai) of Fig. O. A visual inspection of this CMD confirms that 
NGC 1856 harbours an eMSTO. Furthermore, the MS reveals two 
main components below the TO, which are separated by ~0.1 mag 
at raF 336 w ~20.0 and merge together ~1.5 magnitudes below. These 
features are highlighted by the Hess diagram in panel (a 2 ). 

In the CMD of Fig. [3^3 we have arbitrarily selected two sam¬ 
ples of bright MSTO (bMSTO) and faint MSTO (fMSTO) stars 
that we have colored in blue and red, respectively. Moreover we 
have used magenta and green color codes to represent two groups 
of blue MS (bMS) and red MS (rMS) stars. 

To further demonstrate that the broadened MSTO and the dou- 
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Figure 2. Left panels: Original mpssew vs. mpssew ~ ™f 814 W CMD of stars within the NGC 1856 cluster field (green circle in Fig.[T] up), and CMD corrected 
for differential reddening (bottom). The red arrows indicate the reddening direction. Right panels: Original CMD of stars in six sub-regions, I-VI, within the 
cluster field (up). Map of differential reddening for the cluster field. The levels of gray indicate the differential reddening as in the scale plotted on the bottom 
of the figure. The six sub-regions are highlighted with red circles. 


ble MS are not due to differential reddening we have introduced 
the mF438W + ™F814W “ 1.838mF33fjw vs. (mF33(jw “ mF438w)- 1-0065 
(mF438w - niF 555 w) showu in the panel (bi) of Fig. [3 The pseudo¬ 
color and pseudo magnitude used in this diagram are constructed 
in such a way that their total extinction equals to zero, thus making 
the diagram of Fig. [3l)i reddening free, fndeed, by using the ab¬ 
sorption coefficients provided by Aaron Dotter and listed above we 
derive: Af438w + Af814w “ 1-838Af336w=(Af336w - Af 438 w)-1-0065 
(Af 438 w “ Af 555 w)= 0 - The fact that the MSTO and the upper MS 
of NGC 1856 are more broadened than what is expected from ob¬ 
servational errors further demonstrates that the broadening can not 
be due to residual differential reddening. As noticed by the ref¬ 
eree, the Hess diagram and the pseudo-CMD of panels b 2 and b 3 
may reveal more than two sequences. Below mF 438 w + utfshw - 
1.838mF336w~16.2, the rMS and the bMS are clearly visible, while 
at brighter luminosities, around the MSTO, there is some hint of 
a triple sequence. Unfortunately the connection between the MSs 
and the MSTOs is not clear from the diagrams. A visual inspection 


of the Hess diagram of panel b 2 suggests that the rMS could fur¬ 
ther split into distinct MSs, which evolve into the bMSTO and the 
fMSTO. In this case the bMS should evolve into the bluest region 
of the bMSTO. As an alternative the bMS could cross the rMS and 
evolve into the fMSTO. In the next Sect. l5.1l and l5.2l we will further 
investigate these hypothesis. 


3.2 Spatial variations of photometric zero point 

Apart from observational errors, spatial variations of the photo¬ 
metric zero point either due to differential reddening or to small 
uncertainties in the PSF model can broaden the sequences in the 
CMD and mimic multiple stellar populations (e.g. Anderson et 
al. 2008). To investigate this phenomenon, we exploit the mF 438 w 
vs. n 7 F 438 w—^F 8 i 4 w and mF438w vs. nrF 438 w — ^^F 656 N CAfDs and Hess 
diagrams plotted in the lower panels of Fig. [3] The comparison of 
the observed CMD of both panel (ci) and (di) with the error bars 
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™F438W -™F814W ™ F438W ™ F814W 


™ F438W ™ F656N ™ F43BW ™ F656N 


Figure 3. Panel (ai): mpssew vs. mF 336 w ~ ™f 555 W CMD corrected for differential reddening for stars in the cluster field. Panel (a 2 ) shows the corresponding 
Hess diagram for stars around the MSTO. In panel (a 3 ) we have identified a group of bMS and a group of rMS stars that we have colored magenta and green, 
respectively, and a sample of bMSTO and a of fMSTO colored blue and red, respectively. The same color codes have been used to represent the same stars 
in panels (b 3 ), (C 3 ), and (d 3 ). Panels (bi_ 2 , 3 ), (ci,2,3). and (di, 2 , 3 ) are the same as panels (ai_ 2 , 3 ) but for the reddening-free mF438W + tnF 8 l 4 W ~ 1.838mF336w 
vs. (mF336W ~ ™f438w)- 1.0065 (mF438W ~ '«F555w) pseudo CMD, and for the mF438W vs. mF438W “ »tF8i4W and mF438W vs. mF438W ~ »iF656N CMDs corrected 
for differential reddening. 


plotted on the left reveals that the upper MS of NGC 1856 is more 
broadened than expected from photometric errors only. 

To demonstrate that this broadening is not entirely due to pho¬ 
tometric zero-point spatial variations, in the panels (C 3 ) and (d 3 ) of 
Fig.[3we have used the same colors introduced in Fig . [^3 to rep¬ 
resent the bMSTO, fMSTO, bMS, and rMS stars defined therein. 
As suggested by Anderson et al. (2009), if the color and magnitude 
broadening of the MS is entirely due to photometric errors, a star 
that is redder than the bulk of MSTO stars in niF 336 w“”iF 555 w has the 
same probability of being either redder or bluer in mF 438 w “ ™f 814 w 


or mF 438 w - ™f 656 n- Instead we notice that the two groups of bM¬ 
STO and fMSTO stars selected in the CMD of panel (a 3 ) are also 
well separated in the CMDs of panels (C 3 ) and (d 3 ). This further 
demonstrates that the broadening of the MSTO can not be entirely 
due to measurement errors. Similarly, from the mF 336 w - 
the mF 438 w ~ ^F 8 i 4 w> and the mF 438 w “ ^f 656 n colors of the bMS 
and rMS stars we conclude that NGC 1856 hosts at least two MS 
components fainter than the TO. 
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Figure 4. Upper panels: mpssew vs.mF 336 w “ 'Wpsssw CMD of stars in the 
cluster field (left) and for stars in the reference field (right). Lower panels: 
CMD for cluster-field stars after field stars have been statistically subtracted 
(left) and CMD of the subtracted stars (right). 


3.3 Field stars 

Contamination from background and foreground stars can be re¬ 
sponsible for multiple sequences in the CMD. In this subsection 
we present the decontaminated cluster CMD to show that the ob¬ 
served eMSTO and the split MS can not be explained by field stars. 
The best way to derive a CMD solely made of NGC 1856 stars is by 
separating field stars from cluster members using stellar motions. 
Unfortunately, the short baseline of our observations coupled with 
the large distance of the LMC prevents us to derive proper motions 
that are accurate enough to reach this goal. As an alternative, we 
note that, due to the small size of the analyzed WFC3/UVIS field 
of view, we can assume that the distribution of field stars is homo¬ 
geneous. Therefore, to investigate whether the broadened MS of 
NGC 1856 is due to field-star contamination or not, we follow the 
recipe of Paper I and II and compare the CMD of stars in the clus¬ 
ter field and in a reference field which cover the same area as the 
cluster field but includes only field stars. This procedure is illus¬ 
trated in Fig. |4] where we show in the upper panels a zoom of the 
m-pum vs. mF 336 w - mpsssw CMD around the MSTO of NGC 1856 
for stars in the cluster field (upper-left panel) and in the reference 
field (upper-right panel). 

We find that, in the interval of color and magnitude shown in 
Fig.|4]with -0.34 < mF 336 w - mpsssw < 0.47 and 16.7 < mF 336 w < 
21.6 , the cluster and the reference fields contain 2031 and 107 
stars, respectively. This fact suggests that the field contamination 
in the CMD of NGC1856 is of the order of only 5%. 


Then we have statistically subtracted the stars in the reference 
field from the CMD of stars in the cluster field. To do this we have 
defined for each star (i) in the reference field a distance 

~ ^^^FSSSW.cf) — (^F336w rf ^ ^F555W rf hF t^F336W,cf — ^f 336W cf F ’ 

where mF 336 w,cf and raFsssw.cf are the magnitudes of stars in the clus¬ 
ter field, while the magnitudes of stars in the reference field are in¬ 
dicated as mF 336 w,rf and rap 555 w,rf- To account for the fact that the 
color of a star is better constrained than its magnitudes (see discus¬ 
sion in Gallart et al. 2003) we have multiplied the color by a factor 
k=3.2, whose value has been determined as in Marino et al. (2014, 
see their Sect. 3.1). Finally, we flagged the closest cluster-field star 
as a candidate to be subtracted. 

The decontaminated CMD of NGC 1856 plotted in the lower- 
left panel of Fig.|4]demonstrates that the eMSTO and the broad MS 
of NGC 1856 are not due to field stars. For completeness, we plot 
the CMD of subtracted stars in the lower-right panel of Fig.|4] 

3.4 Non-interacting binaries 

Paper I and II have shown that the CMD of all the analyzed young 
and intermediate-age star cluster in the LMC host a large fraction 
of non-interacting binaries, which ranges from ~0.19, in the case 
of NGC 1652 to ~0.46 for NGC 2108. A visual inspection of the 
CMDs of Fig. I^and Fig.j^^i reveals an excess of stars on the red 
side of the MS and suggests that NGC 1856 has a sizeable binary 
population. Since we expect that binaries can affect the distribution 
of colors and magnitudes of MS stars, we have investigated the 
effect of binaries on the morphology of the turn off. To do this we 
first have determined the fraction of binaries along the MS (MS- 
MS binaries). Then we have simulated a CMD of NGC 1856 which 
includes the derived binary fraction and have compared it with our 
observations. 

To infer the fraction of MS-MS binaries in NGC 1856 we have 
used the mF 8 i 4 w vs. mF 438 w “ ™p 814 w CMD, where the binaries with 
q >0.5 are well separated from single MS stars, and followed the 
recipe of Milone et al. (2012). The quantity q is the mass ratio of the 
two components of the binary. Briefly, we have divided the CMD 
into two regions A, and B. Region A includes all the single stars 
and binary systems with a primary star with 20.5< rapguw <21.5, 
and includes both the light- and dark-shaded area of Fig.|3 Region 
B, is the dark-shaded portion of region A and includes binaries with 
q > 0.5. The fraction of binaries with a mass ratio q >0.5 has been 
calculated as 

^>0.5 _ ^CLUSTER ~^REFERENCE _ ^ARTS 

-'BIN _wA 

'^CLUSTER '^REFERENCE '^ARTS 

where A^c(^stfr number of stars, corrected for complete¬ 

ness, in the regions A and B of the CMD of stars in the cluster field 
plotted in the upper-left panel of Fig. [5] ^REPERENCE Corre¬ 

sponding numbers for stars in the CMD of stars in the reference 
field (upper-right panel of Fig. [3, while A^^rts numbers of 

stars in the regions A and B of the CMD made of artificial stars 
(lower-left panel of Fig.|5] see Milone et al. 2012 for details). 

We conclude that NGC 1856 hosts a large fraction of binaries, 
in close analogy with what is observed in other MC clusters. We 
find =0.147±0.015 in the cluster field and infer a total frac¬ 
tion of binaries /g™=0.294±0.031 by assuming a flat mass-ratio 
distribution for binaries with q <0.5. In the lower-right panel of 
Fig.|5]we have used artificial stars and BaSTI isochrone|^ (Pietrin- 
ferni et al. 2004, 2006, 2009) to simulate the mF 336 w vs. rap 336 w “ 

^ http://www.oa-teramo.inaf.it/BASTI 
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^Wpissw CMD of a simple stellar population with age of 320 Myr 
and Z=0.008 (see Sect.O for an age estimate) and by assuming 
the same fraction of binaries above derived. From the comparison 
of the simulated CMD and the observed one plotted in Fig.[^i we 
conclude that non-interacting binaries alone are not able to repro¬ 
duce neither the eMSTO nor the double MS observed in NGC 1856. 

The evidence of an eMSTO in NGC 1856 is in disagreement 
with the conclusions by Bastian & Silva Villa (2013) who have an¬ 
alyzed WFPC2 photometry in F555W and F814W of NGC 1856 
and found no evidence for an eMSTO in this cluster. Their work 
is based on the catalog by Brocato et al. (2001) who provided high 
quality photometry for the standard available at that time. Our find¬ 
ing of an eMSTO and a split MS in NGC 1856 is likely due, in 
part, to the high quality of the photometry obtained from images 
taken with the new WFC3AJV1S and ACS/WFC cameras of HST 
and by using the recent techniques developed by Jay Anderson and 
collaborators (see Sect.j^for details.) 

Moreover, we have shown that the eMSTO is better visi¬ 
ble when we use colors involving the F336W filter (like F336W- 
F555W and F336W-F814W colors colors, see Fig.|2]and|3]l- There¬ 
fore the F336W photometry used in this paper is more efficient than 
F555W and F814W to detect the eMSTO in young clusters. 


4 THE DOUBLE MS OF NGC 1856 

Section [3] has provided evidence that, below the MSTO, the MS 
of NGC 1856 is split, and that such splitting is real. To our knowl¬ 
edge, after the case of NGC 1844 (see Paper 11), this is the second 
time that a double or broad MS in a young star cluster has been 
detected. This finding deserves further analysis to determine the 
fraction of rMS and bMS stars. In order to do this, we have ap¬ 
plied to NGC 1856 a procedure which has been widely used in pre¬ 
vious studies of multiple stellar populations in GCs (e.g. Piotto et 
al. 2007) and is illustrated in Fig.|^ 

The double MS of NGC 1856 is clearly visible in the panel a 
of Fig. 1^ where we show the /npsssw vs. mp 33 jw “ ™f 555 w CMD of 
stars in the cluster field. Only stars within the gray rectangle over 
the magnitude interval 19.9< mpsssw <20.6, where the MS split is 
most prominent, are used in the following analysis. This region of 
the CMD is zoomed in panel b where we also plot the fiducial line 
of the red MS. To determine the line we selected a sample of red 
MS stars by eye and calculated the median color and magnitudes 
of these stars in intervals of 0.2 mag in the F555W band. These 
median points have been interpolated with a cubic spline. The ver- 
ticalized mpjsjw vs. A(mp 336 w ~ ™f 555 w) CMD shown in panel c of 
Fig-|6]is obtained by subtracting the color of the fiducial at the cor¬ 
responding mpjssw from the color of each star. The color separation 
between bMS and rMS increases with luminosity and ranges from 
-0.1 at mp555w=20.6 to ~ -0.15 at mF555w=19.9. The distribu¬ 
tion of A(mp 335 w - J?!p 555 w) is plotted in panels d in four bins of 
™F 555 W- The histogram distributions exhibit two peaks that we have 
fitted with two Gaussians (magenta for the bMS and aqua for the 
rMS). From the area under the four Gaussian profiles we infer that 
the bMS and the rMS host 33±5% and 67±5% of MS stars, respec¬ 
tively, in the selected magnitude interval. 


5 POSSIBLE ORIGIN OF THE MS SPLIT AND OF THE 
EMSTO 

In order to shed light on the physical reasons at the basis of the 
double MS and the eMSTO of NGC 1856, we have compared the 
observed CMD with isochrones from the BaSTI archive. To do this, 
we have used isochrones which account for the core-convective 
overshoot during the central H-burning stage as in Pietrinfemi et 
al. (2004). Specifically, in Sect. 15. l| we test the hypothesis that the 
spread around the MSTO and in the MS are due to age variation 
only, while, in Sect. l5.2l we investigate possible etfects of C-fN-i- 0 
and iron variation on the CMD morphology. 

5.1 Age and age spread 

Several authors suggested that the eMSTO of intermediate-age 
clusters in the MCs is due to a prolonged star-formation history 
(e.g. Mackey et al.2008; Paper I; Goudfrooij et al.2011). In this 
section we investigate the hypothesis that the eMSTO of the young 
cluster NGC 1856 is due to multiple stellar generations with the 
same chemical composition and infer the consequent spread in age. 

Figure |3 shows that, adopting a distance modulus (m - 
M)o=18.35, a reddening £(6-V)=0.15, a metallicity Z=0.008, and 
[a/Fe]=0 the average distribution of stars in the mp 336 w vs. mF 336 w- 
™f 555 w and in the raF 438 w vs. rap 438 w “ ™f 814 w CMDs can be repro¬ 
duced by assuming a 300 Myr age, while the blue and the red en¬ 
velopes of the MS and the MSTO are well matched by a 250-Myr 
and 400-Myr old isochrone, respectively. 

The comparison between isochrones and observations pro¬ 
vides a rough estimate of the duration of the star formation in 
NGC 1856 (in the hypothesis of prolonged or multiple star for¬ 
mation). To determine a more accurate age distribution of stars in 
NGC 1856 we have used the procedure illustrated in Fig. [ 8 ] We 
started by superimposing two red lines on the mF 336 w vs. mp 33 jw - 
'npsssw CMD, as shown in panel (ai). These two lines have been 
drawn by hand with the criterion of selecting the region around 
the MSTO where the color and magnitude spread due to age varia¬ 
tion is more clearly distinguishable. Only the stars between the two 
lines and belonging to the decontaminated CMD of Sect. l3.3l (black 
circles in Fig.[ 8 jii) have been used to estimate the age distribution. 

To do this, we have superimposed on the CMD the grid of 
isochrones available from the BaSTI database for the same values 
of metallicity, [a/Fe], distance and reddening as above described. 
The available isochrones from the BaSTI webpage - represented in 
Fig.[8]with thick gray lines - have ages between 100 Myr and 550 
Myr in steps of 50 Myrs. The isochrones represented with thin gray 
lines are separated by 10 Myr and have been obtained by linearly 
interpolating the BaSTI isochrones. 

In order to estimate the age of a given star, j, we started 
by calculating the distance in the CMD between star j and all 
isochrones. Then we have defined the quantities D, =d|/(di -l-d 2 ) and 
D 2 =d 2 /(di-l-d 2 ), where di and d 2 are the distances from the two 
closest isochrones of age ti and t 2 of a given star. The resulting best 
age estimate is: fj = 1^02 + t 2 Di 

To derive the kernel-density distribution of ages for the ana¬ 
lyzed stars, we assumed a Gaussian kernel with dispersion equal 
to the observational error a. Then we have divided the age interval 
100<t<550 Myrs into bins of 1 Myr and, for each bin, we have cal¬ 
culated the quantity ipi = j-exp ^ where cj is the com¬ 

pleteness of the star), and NTOT=1914 is the number of analyzed 
stars between the two red lines. The normalized age distribution is 
represented with a black line in panel (a 2 ) of Fig.[^ 



8 A. P. Milone et al. 



Figure 5. mpgi4w vs.mF438w - mpSMW CMD for stars in the cluster field (upper-left), in the reference field (upper-right), and simulated CMD (lower-left). 
Light- and dark-gray areas mark the region A of the CMD used to estimate the binary fraction. Region B is the sub-region of A colored dark-gray. Candidate 
single stars and binaries with q<0.5 in the analyzed interval of magnitude are represented with black points, while crosses represent candidate binaries with 
q>0.5. The simulated mp 335 w vs. mp 33 (jw “ 'ttpsssw CMD is plotted in the bottom-right panel. See text for details. 


Each isochrone of age t intersects the two red lines of Fig.at 
two points, Pb,t and Pf t, which correspond to mass Mh and Mf. As 
an example, in Fig.[^ we have indicated with large red circles the 
points corresponding to the 200-Myr and 450-Myr old isochrones. 
The mass interval dA1=Alb - Mf changes from one isochrone to 
another as shown in the panel (a 3 ) of Fig. [8] where we plot AM as 
a function of the age of the isochrone. 

As a consequence, the distribution 0 ; above obtained clearly 
depends on our choice of the CMD region used to determine it. 
To remove such dependence from the derived age distribution, we 

have calculated the quantity c AM- ’ "'here 

dM^ is the mass interval corresponding to star j. The normalized 
age distribution is represented with a red line in the panel (a2) of 
Fig.i 

In order to estimate the observed age spread we have calcu¬ 
lated the quantities ageio and agecio defined as the age values that 
satisfy the conditions: 2 -floo, <^* = 2 f=foo, ^,* and loo, 0 * = 
M 2 ?=i’oo, respectively, and defined the quantity o-8o,oBS=age9o- 
ageifl. Similarly, we have defined the quantity (T 68 ,obs- We obtained 
o' 8 o,obs = 185 Myr and fr 68 ,OBS = 146 Myr. 

Note that the age spread above estimated is not the intrinsic 
age dispersion. In fact the measured age spread (cr80(68),OBs) results 
from the combination of the intrinsic age dispersion ((T80(68),age) 


and the broadening due to measurement errors and binary stars 

(rr80(68 ),err); f^80(68),age = 

To infer (T 80 ( 68 ),err and determine the intrinsic age spread we 
have used AS tests and simulated the mp 336 w vs. mp 336 w - mpsssw 
CMD plotted in panel (bj) of Fig. [8] We have followed the recipe 
described in detail in Paper I (see their Sect. 6.2). We used the same 
metallicity, a abundance, distance modulus and reddening adopted 
for the isochrone fit of Fig. |7] and the same fraction of binaries as 
derived in Sect. 13.41 We assumed the age distribution represented 
by the black histogram in panel (bg) of Fig. [8] (the dashed black 
line represents the corresponding kernel-density distribution of the 
input ages for our simulation). 

We have applied to the simulated CMD the same procedure as 
used for the observed CMD to determine the age spread and derived 
the red line of panel (bg). Note that the age distribution adopted 
as input to generate the simulated CMD [black histogram in panel 
(bg)] has been derived with the criteria that the kernel-density age 
distribution derived from the simulated diagram [red line in panel 
(bg)] reproduces the age distribution obtained from the observed 
CMD [(ag), red line]. We emphasize that the detailed morphology 
of the histogram representing the input age distribution for the sim¬ 
ulated CMD could have been different. Here, we are interested in 
estimating the intrinsic age dispersion (T68,age and crgo.age not the 
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Figure 6. mpssjw vs. mpssew ~ w*F 555 W CMD of stars in the cluster field (panel a). Panel b shows a zoom of the region where the MS split is more clearly 
visible. The red line overimposed on the CMD is the fiducial line of the rMS. The verticalized mpsssw vs. A(mF 336 w “ wtF 555 w) diagram is plotted in panel c, 
while panels d show the histogram of the A(mp 335 w ~ 'ttF 555 w) distribution for stars in the four magnitude intervals delimited hy gray-horizontal lines in panels 
b and c. The least-squares best-fit bi-Gaussian function is colored gray, while the two components are represented with magenta and aqua color codes. 




■™P336W ™F555W F438W ™F814W 

Figure 7. The three isochrones from the BaSTI database which provide the best match with the observed CMDs. Specifically, the blue and the red isochrones 
match the bluer and redder envelope of the eMSTO of NGC 1856, respectively, while the green ones fit the average distribution of cluster stars in the CMD. 
See text for details. 


detailed distribution of the ages of the stars in NGC 1856. What 
is relevant from our analysis is the similarity of the observed and 
simulated red lines in panels (a 2 ) and (b 2 ) of Fig.[8] 

Also for the simulated CMD we defined a ctso.sim and crsg sim^ 
in close analogy with what was done for the observed CMD. We 
obtain cr8o,siM=182 Myr and cr68,siM=145 Myr. Similarly, the dif¬ 
ference between the 90* and the 10* and the 84* and the 16* per¬ 
centile of the age distribution used as input of our simulation [black 
histogram of panel (b 2 )] have been called cr8o,iNP and cr68,iNP, re¬ 
spectively. We have found cr8o,iNP=167 Myr and (T58,inp=138 Myr. 
Therefore the contribute of measurement errors and binary stars 
to the observed age spread is <rso,ERR = t/'^so.sim -^o.w^'^ 2 Myr and 
d'68.err = “ ‘T8.inp“ 45 Myr. 

At this point we have estimated the intrinsic age dis¬ 


persion as T8„,.ge = V'^go,oBs-‘^8o,Effi=l'74±31 Myr and = 

V°'68 0 Bs “‘T8 err“ 139±22 Myr. Uncertainties on cr68 and ctso are esti¬ 
mated by bootstrapping with replacements performed 1,000 times 
on both the observed and the simulated age distributions. The er¬ 
rors indicate one standard deviation (68.27"’ percentile) of the boot¬ 
strapped measurements. 

5.2 C-i-N+O and metallicity. 

Multiple stellar populations in nearly all Galactic GCs are char¬ 
acterized by star-to-star variation in several light elements, like 
C, N, Na, and O (e.g. Kraft 1994; Gratton et al. 2004, 2012; Car- 
retta et al. 2009). Moreover stars with different light-element abun¬ 
dance populate distinct sequences in CMDs built with ultraviolet 
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Figure 8. This figure illustrates the method to derive the age spread of NGC 1856. Panel (ai) reproduces the mpjsew vs. mpssew — ™f 555 W CMD of Fig.|2] 
with overimposed a grid of isochrones from the BaSTI database. Isochrones with ages from 200 to 500 Myrs in steps of 50 Myrs are represented with thick 
gray lines, while the isochrones plotted with thin lines are separated hy 10 Myrs. The two red lines delimit the CMD region where the color and magnitude 
broadening is more evident. Only stars in this region, which are plotted with small black circles, are used to derive the distribution of ages. We have indicated 
with large red circles the intersections (Pb,f, 200 . and Pb,f, 45 o) between the 200- and 450-Myr old isochrones and the two lines. The kernel-density distribution 
of the ages plotted in panel (a 2 ) and derived for stars between the two red lines is represented with a black line, while the red lines represent the global 
distribution inferred for NGC 1856. In panel (as), we have plotted the mass difference dAI as a function of the age (see text for details). Panel (bi) shows the 
simulated mpssew vs. mpssew ~ ™F 555 W CMD obtained by assuming the age distribution corresponding to the black histogram of panel (bs). The black-dashed 
line indicates the kernel-density distribution of the ages of simulated stars, while the red-continuous line is the distribution inferred from the analysis of the 
simulated CMD of panel (bi). 


magnitudes (e.g. Marino et al. 2008; Yong et al. 2008; Sbordone et 
al. 2011; Dotter et al. 2015). In some massive Galactic GCs, there 
are also internal variations in iron and in the overall C-fN-l-0 abun¬ 
dance (e.g. Norris & Da Costa 1995; Cassisi et al. 2008; Marino et 
al. 2009; Ventura et al.,2009; Yong et al. 2014) which are respon¬ 
sible for multiple sequences in the CMD (e.g.Piotto et al. 2005; 
Marino et al. 2011; Milone et al. 2015). 

In absence of a detailed study on the chemical-composition 
of NGC 1856, in this subsection we investigate a possible effect of 
C-l-N-l-O and iron variation on the CMD morphology. First of all 
we note that a visual inspection of the mF 555 w vs. mF 33 (jw “ ™f 555 w 
H ess diagram of stars in the cluster field plotted in the left panel of 


Fig. |9] could suggest that the rMS and bMSTO are connected and 
that the bMS evolve into the fMSTO. 

The middle panel of Fig.l^shows that a solar-scaled isochrone 
for Z=0.01, [Fe/H]=-0.25, and age, t=300Myr, matches the bulk 
of stars in NGC 1856 and reproduces the rMS and the bMSTO. The 
same figure shows that the bMS and the fMSTO are consistent with 
solar-scaled, metal poor isochrones with Z=0.005, [Fe/H]=-0.6 
and age of 400-450 Myr, although the fit is quite poor. 

In the right panel of Fig.j^we have fitted the rMS and the bM¬ 
STO with an nr-enhanced isochrone with Z=0.015, [Fe/H]=-0.7, 
age of 275 Myr and with a factor of ~2 enhancement in CNO- 
element sum (blue line). The red isochrones correspond to a- 
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enhanced models for Z=0.005, [Fe/H]=-0.6 and age of 400- 
450 Myr. 

In summary, the comparison between the observed CMD and 
the isochrones of Fig. [9] shows that the bMS and the fMSTO are 
well fitted by an old population with [Fe/H]=-0.6, solar-scaled 
CNO, and age of ~400Myr, while the rMS and the bMSTO are 
consistent with a ~300 Myr-old stellar population. The later can be 
either enhanced in [Fe/H] by ~0.4 dex or can have almost the same 
iron content as the old population but is enhanced in C-l-N-l-0 by 
about a factor of two. In all cases, even assuming variations in C, 
N, O and Fe we need an age spread of ~ 150 Myr in order to repro¬ 
duce the broadened distribution of the upper MS and eMSTO. 

Finally, we enphasize that the analysis above only provides an 
attempt to interpret the eMSTO and the split MS of NGC 1856 and 
that the present dataset does not allow us to to establish whether 
NGC 1856 has homogeneous C-l-N-rO and iron. Spectroscopy of 
NGC 1856 stars is mandatory to properly contraint the chemical 
composition of this cluster. 

5.3 The extended red clump. 

Left panel of Fig. [li shows a zoom of the mF 336 w vs. mF 336 w - 
"tFassw around the red clump. The 142 black points represent stars 
in the cluster field, while the 8 reference-field stars are plotted with 
red crosses. The contamination from field stars is lower than ~6%. 
We have used the same procedure of Sect. l3.3l to statistically sub¬ 
tract field stars from the cluster-field CMD and plot the CMD of 
subtracted stars in the middle panel of Fig. The red clump is 
widely spread in color and magnitude and spans a range of ~0.4 
mag in mF 336 w - "Jfsssw and show same hint of bimodality, with 
two main components at (raF 336 w - raF555w:'UF336w)=(L50:19.25) 
and (1.60:19.25). We have used solar-scaled BaSTI isochrones with 
Z=0.01, and [Fe/H]=-0.25, to simulate a simple stellar population 
with t=300Myr. We assumed the same errors as in the observed 
CMD, and the fraction of binaries measured in Sect. l3.4l The simu¬ 
lated CMD is plotted in the right panel of Fig. where we plotted 
the same number of stars as in the middle-panel CMD. A visual 
comparison of the observed and the simulated CMD immediately 
reveals that the red clump of NGC 1856 is not consistent with a 
simple population. 

In principle, the red clump could be used to constrain the star- 
formation history of NGC 1856. However, the luminosity and the 
color of the red clump depends not only on age but also, on several 
parameters like the amount of core overshooting, the mass loss, 
and the helium content (see e.g. Girardi et al. 2009). Unfortunately, 
these quantities are not constrained for NGC 1856, and we prefer 
to not use the red clump to investigate the star-formation history of 
NGC 1856. 


6 SUMMARY AND DISCUSSION 

We have used multi-band UVIS/WFC3 photometry to study the 
LMC star cluster NGC 1856. We found that the MSTO is broad¬ 
ened, and that, below the turn-off, the MS splits into two main 
components. In the magnitude interval 19.9< ujfsssw <20.6 the 
blue and the red MSs host 33±5% and 67±5% of the total number 
of MS stars, respectively. 

We have demonstrated that both the split MS and the eMSTO 
can not be due to photometric errors, differential reddening, or un¬ 
resolved binaries. As discussed in Sect.[T] the eMSTO is a common 
feature of intermediate-age star clusters in the Large and the Small 


Magellanic Cloud although its interpretation is still controversial. 
This paper shows the first evidence of an eMSTO in a young star 
cluster with age less than ~1 Gyr and provides new steps towards 
the understanding of this phenomenon. We have compared the ob¬ 
served CMD with isochrones from the BaSTI archive and estimated 
an average age of ~300 Myr. We have demonstrated that, if stars in 
NGC 1856 are chemically homogeneous and the eMSTO is due to 
a prolonged star formation, the latter should have lasted ~150 Myr, 
not necessarily in a continuous way, in view of the features we ob¬ 
serve in the upper MS and MSTO. Specifically, the interval between 
the 90* and the 10* percentile of the age distribution corresponds 
to CTgo = 174±31 Myr. The age difference between the 84* and the 
16* percentile is ixeg = 139±22 Myr. In the hypothesis that the MS 
and MSTO broadening are due to a prolonged star formation, the 
red MS evolves into both the faint and the bright MSTO, and the 
blue MS into the bright MSTO. 

The eMSTO of intermediate-age MC clusters have been inter¬ 
preted either with an age spread or with stellar rotation. The idea 
that this feature of the CMD is due to prolonged star formation has 
been suggested by several authors (e.g. Mackey et al. 2008; Glatt et 
al. 2008; Paper I; Goudfrooij et al. 2011) and has been recently sup¬ 
ported by Goudfrooij et al. (2014) on the basis of their study of 18 
MC clusters. These authors have found a correlation between the 
MSTO width and the escape velocity (Vesc) of the host cluster and 
argued that these results are unlikely due to stellar rotation. 

Goudfrooij et al. (2014) and Correnti et al. (2014) noticed that 
all the analyzed clusters have Vesc > 17 km s“' and suggested 
that eMSTO can occur only in clusters whose escape velocities 
are higher than the winds of first-generation stars from which the 
second generations formed. Intermediate-mass asymptotic-giant- 
branch stars and massive binaries have wind velocities close to 12- 
15 km s“', thus being possible candidates for producing the ma¬ 
terial from which the second generation formed (see Goudfrooij 
et al. 2014 and Correnti et al. 2014 for details). Moreover, Goud¬ 
frooij, Correnti and collaborators have hypothesized the existence 
of an early-escape-velocity threshold of about 15 km s“* which dif¬ 
ferentiate clusters that have experienced a prolonged star-formation 
from those without eMSTO. 

The evidence of an eMSTO in NGC 1856, whose velocity es¬ 
cape did not exceed ~15km s“' at its formation (Goudfrooij et 
al. 2014), is consistent with this picture. In this context, it would 
be interesting to understand why the ■''180-Myr old LMC clus¬ 
ter NGC 1866, whose velocity escape is only slightly smaller than 
that of NGC 1856 does not show an eMSTO (Bastian & Silva- 
Villa2013). 

We have investigated the possibility that the double MS and 
the eMSTO have different metallicity and C-rN-t-O abundance. In 
these cases the blue MS would evolve into the faint MSTO, while 
bright-MSTO stars would be the progeny of the red MS. How¬ 
ever, also in case of iron variation among stars in NGC 1856, we 
need to assume an age difference of ~150 Myr. If we assume 
that NGC 1856 has constant overall CNO content, the blue MS 
and the faint MSTO correspond to a ~400Myr population with 
[Fe/H]=-0.6, while the rMS and the bMSTO are consistent with 
a ~150Myr-younger stellar population enhanced in [Fe/H] by ~0.4 
dex. Observations are also consistent with the younger population 
having almost the same iron content as the old ones but being en¬ 
hanced in C-t-N+O by about a factor of two. 

Interestingly, in Paper II we have found that the ~ 150-old 
LMC cluster NGC 1844 exhibits an intrinsic color spread along the 
MS. We have shown that such a feature is well fitted by isochrones 
with different C-i-N-hO, although these models provide a poor fit of 
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Figure 9. Leji panel: Hess diagram of the mpsssw vs. wf336W “ ^F555 W CMD of Fig.[^ Middle panel: Solar scaled isochrones with different metallicity and 
age from the BaSTI database are superimposed on the observed CMD. Right panel: As in the middle panel but for alpha-enhanced isochrones. In this case the 
isochrone with Z=0.015 is also CNO enhanced. See text for details. 
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Figure 10. Left panel: mp 335 w ''s. mp 335 w - mpsssw CMD zoomed around the red clump. Black circles and red crosses are stars in the cluster and reference 
field, respectively. Middle panel: CMD of stars in the cluster field after that reference-field stars have been statistically subtracted. Right panel: Simulated 
CMD of a 300 Myr-old simple stellar population (see text for details). 

the bright MS of NGC 1844. Due to the possible similarity between 
the MS of NGC 1844 and NGC 1856 it is tempting to speculate that 
C-l-N-t-0 variations could play a role on the morphology of the MS 
and MSTO morphology of young and intermediate-age MC clus¬ 
ters. 

Unfortunately no abundance measurements are available for 
NGC 1856 to date. Moreover, we are not able to establish whether 
NGC 1856 has homogeneous C, N, O and Fe from photometry 
alone because the connection between the faint and the bright 
MSTO with the two MSs can not be established unequivocally 
from our CMDs. In this context, it is worth noting the finding 
by Mucciarelli et al. (2014) that in the eMSTO of the LMC clus¬ 
ter NGC 1846 there is no evidence for significant variation of 
[O/Fe], [Na/Fe], [Al/Fe], [Mg/Fe], and [Fe/H]. This result suggests 


that intermediate-age clusters in MCs could have homogeneous 
iron and light-element abundance. A chemical-abundance study of 
NGC 1856 is needed in order to estimate the iron and C, N, O con¬ 
tent of this cluster. 


As an alternative interpretation of the eMSTO phenomenon, 
Bastian & De Mink (2009) have suggested that rotation in stars 
mimic an age spread and can be responsible for the eMSTO of 
intermediate-age MC clusters. This idea is supported by the fact 
that the structure of a star is significantly affected by rotation and 
that its effective temperature depends on the inclination angle rel¬ 
ative to the observer. This interpretation however is in disagree¬ 
ment with the results by Girardi et al. (2011), who have calculated 
isochrones for intermediate-age MC clusters and concluded that a 
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dispersion of radial velocity can not be the only reason for the eM- 
STO. 

Platais et al. (2012) have measured projected rotational veloc¬ 
ities for upper-MS stars of the ~1.3 Gyr old Galactic open clus¬ 
ter Trumpler20 and found very-high rotations up to v sini~180 
km s-1. These authors have selected two groups of slow- and fast- 
rotators and found that fast rotators have a marginal bluer colors 

than slow rotators (d(V - /)-0.01 mag). They concluded that the 

MSTO morphology of Trupler 20 is marginally affected by rota¬ 
tion. Unfortunately no direct measurements of rotational velocities 
are available for MC clusters with eMSTO. 

An accurate comparison of our observations with appropriate 
stellar models is certainly needed to understand the importance of 
stellar rotation for the explanation of the eMSTO in NGC 1856. 
Nevertheless our findings provide a major challenge to the scenario 
suggested by Bastian & De Mink. Indeed, in order to explain the 
eMSTO of intermediate-age clusters, these authors have suggested 
that only stars with masses \.2 < M < I.IMq would be affected 
by rotation. In the young NGC 1856, stars in this mass interval are 
located on the MS region, with 20.7< mpssew <23.4, where there is 
no evidence for a large spread in color, thus suggesting that stellar 
rotation, if present, would not significantly affect the colors and 
the magnitudes of these stars. If the eMSTO of both NGC 1856 
and intermediate-age clusters is due to stellar rotation it would be 
challenging to understand why this phenomenon would affect stars 
with different masses in clusters with different ages. 

Yang et al. (2011) have suggested merged binary systems and 
interactive binaries with mass transfer can be responsible of both 
the eMSTO and dual clumps in intermediate-age MC clusters. De¬ 
spite this the rest of the MS of these clusters is not significantly 
spread by binary interaction. The evidence that NGC 1856 exhibits 
a dual MS seems in disagreement with the scenario suggested by 
Yang and collaborators. However, it should be emphasized that con¬ 
clusions from these authors are based on intermediate-age star clus¬ 
ters and that appropriate studies of young clusters are mandatory to 
establish if interacting binaries are consistent with observations of 
NGC 1856 or not. 

In summary the detection of a double MS and an eMSTO 
in the young NGC 1856 adds new information to the puzzling 
eMSTO-phenomenon of MC star clusters. On the one hand, the 
new findings seem consistent with a prolonged star-formation his¬ 
tory thus challenging the interpretation that either stellar rotation or 
interacting binaries can be responsible of the eMSTO. On the other 
hand, the recent findings that the SGB morphology of some massive 
intermediate-age MC cluster with eMSTO is not consistent with a 
spread in age (Bastian & Niederhofer2015; Li et al. 2014), and the 
lack of evidence of prolonged star-formation history in extragalac- 
tic massive clusters raise doubts on the interpretation of eMSTO as 
due to multiple stellar generations with different ages (or prolonged 
star formation). 
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